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Abstract. We have analyzed the emission line properties 
in the integrated spectra of 15 normal spiral galaxies. We 
show that very clear trends appear when plotting rele- 
vant emission line ratios or equivalent widths as a func- 
tion of galaxy spectral types, obtained with a Principal 
Component Analysis of the continua and absorption fea- 
tures of spectra. The equivalent widths of all the lines 
analyzed correlate extremely well with spectral types, im- 
plying that each of them can be considered a good indica- 
tor of the spectral type in normal galaxies. The position 
of most galaxies of our sample in classical emission line 
diagnostic diagrams follows that of individual giant HII 
regions in spiral galaxies, but for the earliest type galaxies, 
the emission line pattern resembles more that of LINERs. 
Therefore, the direct interpretation of equivalent widths 
in terms of star formation rates would be misleading in 
such cases. The observed trends in the emission line ratios 
as a function of galaxy spectral type suggest a decrease 
of O/H, a decrease of N/O, an increase of the average 
effective temperature or ionization parameter, and a de- 
crease of the effective internal extinction of galaxies with 
increasing (early to late) spectral type. 

Key words: methods: statistical - galaxies: fundamental 
parameters - galaxies: ISM - galaxies: abundances - galax- 
ies: spiral 



1. Introduction 

Emission lines in normal galaxies are powerful tracers of 
the physical processes associated with the galaxy inter- 
stellar medium and the formation of massive stars. They 
allow quantitative estimates of the star formation rates, as 
well as of the chemical abundances and physical conditions 



of the gas in galaxies (McCall et al. 1985; Gallagher et al 



1989; Kennicutt, Tamblyn & Congdon 1994 Zaritsky et 



al. 1994; Barbara & Poggianti 1997) 



A census of emission line properties in integrated spec- 
tra of normal galaxies is particularly interesting in view 
of studies of galaxies at large redshifts, too distant to be 
spatially resolved. It provides useful information also for 
nearby galaxies, as it allows investigating the global prop- 
erties of their emission. It is well known that emission lines 
tend to become more prominent as one goes from early to 
late galaxy types, and the degree to which the spiral arms 
are resolved into individual HII regions (where most of the 
emission lines in normal galaxies are produced) is one of 
the most important criteria for discriminating among the 
spiral types in the Hubble system (Hubble 1936; Sandage 
1961). Kennicutt (1992a, hereafter K92a) has already dis- 
cussed some properties of the emission lines in integrated 
spectra, mainly evaluating their reliability as quantitative 
tracers of the total massive star formation rate in galaxies. 



Lehnert & Heckman (1994) have examined the location of 
integrated spectra of galaxies in the standard emission line 
diagnostic diagrams used to classify emission line objects 
(Bald win, P hillips & Terlevich |1981| ; Veilleux & Oster- 
brock 1987 ), concluding that the integrated emission in 
star-forming galaxies has a substantial contribution of a 
diffuse component, with physical properties different from 
those found in high surface brightness HII regions. Most 
studies, however, have focused on the properties of emis- 
sion lines produced either in HII regions in the spiral arms 
of disk galaxies (e.g. M cCall et al. 1985 ; Be lley fc Roy 



1992 



Oey & Ke nnicut t |1993| ; Zarits ky et al. [1994t Ken- 
nicutt & Garnett |1996| ; Roy fc Wal sh |l997| ) or in g alacti c 
nuclei (Heckman et al. |l980| Keel |1983| ; Ho et al. |1997| ). 
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In particular, Zaritsky et al. ( 1994 ) have investigated the 
oxygen abundance properties of a sample of disk galax- 
ies from the spectra of individual HII regions located at 
various galactocentric radii. They found that the charac- 
teristic abundance of the galaxies correlates well with both 
their morphological type and their luminosity. 

In this paper we examine the trends of emission 
line properties in the integrated spectra of 15 normal, 
nearby spiral gala xies. T he spectra come from Kennicutt 's 
(1992b, hereafter K92b ) spectro-photometric atlas. 
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The galaxies are ordered by galaxy spectral type (here- 
after ST), obtained using Principal Component Analysis 
of the integrated properties of the continua and absorp- 
tion features of 23 normal galaxies of all morphologies, in 
a manner similar to the one introduced by Sodre & Cuevas 

Connolly et 
Galaz 



al. 1991), their spectral types (see below), and the spec- 
tral classification of their nuclei according to Ho ct al. 



(1994, 1997) and several other authors (e.g.. 
al. |1995| ; Zaritsk y et a l. |1995| ; Folkes et al. 
& de Lapparent 1998). Such a procedure allows to de- 



1996 



fine a spectral classification that correlates well with Hub- 
ble morphological types and present some advantages over 
the usual morphological classification. Firstly, it provides 
quantitative, continuous, and well defined types, avoiding 
the ambiguities of the intrinsically more qualitative and 
subjective morphological classification. Secondly, the se- 
quence of galaxies obtained in this way is easier to model 
than the Hubble sequence, because the physical process 
behind galaxy spectra are better understood than those 
needed to explain galaxy morphologies. Thirdly, this clas- 
sification can be applied to redshift surveys where no infor- 
mation is available on the galaxy morphologies, allowing 
to "recycle" data obtained for other purposes (de Lappar- 
ent, Galaz & Arnouts |l998t Bromley et al. |1998| ). 

We will show that the patterns of variation of emis- 
sion line properties along the sequence of normal galaxies 
are much more regular when considering the galaxy spec- 
tral type sequence than when using the morphological se- 
quence. 

This paper is organized as follows. In Sect. 2 we in- 
troduce the sample of normal galaxies, briefly describe 
the procedure used to obtain spectral types from inte- 
grated galaxy spectra, and present the relevant data for 
the emission lines in the blue/ visible region of the spectra. 
In Sect. 3 we present the trends of equivalent widths and 
emission line ratios with galaxy spectral types. In Sect. 4 
we interpret these trends, after a brief reminder on emis- 
sion line theory and a comparison of the loci of integrated 
galaxy spectra with those of giant HII regions in classical 
diagnostic diagrams. Finally, in Sect. 5 we summarize our 
results. 



2. The data base 

2.1. The sample of galaxies 

The integrated spectra of galaxies discussed in this work 
come from the atlas of Kennicutt ( K92bj ) . They are suit- 
able for setting the basis of a global characterization of 
galaxies because this is the largest sample of integrated 
spectra of low-redshift galaxies available and obtained in 
a uniform way. From this atlas we, conservatively, rejected 
all the galaxies showing any evidence of strong spectral pe- 
culiarities or interactions (i.e. AGNs, starbursts and merg- 
ers). This left us with a total of 23 normal galaxies, which 
will be our standards. They are presented in Table 1, to- 
gether with their Hubble types (taken from K92a ), their 
T-types (following the RC3 convention- de Vaucouleurs et 



(1997; see Sect. 4). These galaxies are all at high Galactic 
latitude, so foreground absorption by our Galaxy is unim- 
portant. Note that the set of standards adopted here is 
very similar to those used to characterize normal galaxies 
in other studies of spectral classification (e.g. Folkes et al. 



1996; Galaz & de Lapparent 1998) 



2.2. The spectral classification 

Our general approach to spect ral cla ssification is described 
i n det ail in Sodre & Cuevas ( 1997 ) and in Cuevas et al. 
(1998). Here we only present an overview of the proce- 
dures used to obtain spectral types for the galaxies dis- 
cussed here. The point to be stressed is that the spectra 
of normal galaxies form a sequence- the spectral sequence- 
in the spectral space spanned by the M-dimensional vec- 
tors that contain the spectra, each vector being the flux 
of a galaxy (or a scaled version of it) sampled at M wave- 



lengths (Sodre fc Cuevas 1994, 1997; Connolly et al. 1995 



Folkes et al. 1996). The spectral sequence correlates well 



with the Hubble morphological sequence, and we define 
the spectral type of a galaxy from its position along the 
spectral sequence. 

The spectral sequence is identified using Principal 
Component Analysis (PC A). PC A allows to define a new 
orthonormal reference system in spectral space, centered 
at the baricenter of galaxy spectra and with basis- vectors 
(the principal components) spanning directions of maxi- 
mum variance. The plane defined by the first two prin- 
cipal components- which we call principal plane- is the 
plane that contains the maximum variance in the spectral 
space. 

PCA was applied to a pre-processed version of the 23 
original galaxy spectra. Firstly, the spectra were shifted 
to the rest frame and re-sampled in the wavelength inter- 
val from 3784 A to 6500 A, in equal-width bins of 2 A. 
Secondly, we removed from the analysis 4 regions of ~20 
A each centered at the wavelengths of He II A4686, H/3 
A4861, [O III] A4959 and [O III] A5007. This was done in 
order to avoid the inclusion of emission lines in the analy- 
sis, that increases the dispersion of the spectra in the prin- 
cipal plane (mainly due to an increase in the second princi- 
pal component). The spectra, now sampled at M = 1277 
wavelength intervals, were then normalized to the same 
mean flux (%2\ A — !)■ After, we computed a mean spec- 
trum and, finally, we subtracted this mean spectrum from 
the spectrum of each galaxy. PCA was then used to obtain 
the principal components. This procedure is equivalent to 
the conventional PCA on the covariance matrix, that is, 
the basis vectors are the eigenvectors of a covariance ma- 
trix. Note that the principal components are dependent 
of the data and its pre-processing, and using a different 
normalization for the spectra or taking some other wave- 
length interval certainly would lead to a different result 
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Table 1. Standard galaxies and classification 



Name 


Hubble type 


T-type 


spectral type 


nuclear types 


NGC 3379 


E0 


-5 


-5.2 




NGC4472 


E1/S0 


-4 


-5.7 




NGC 4648 


E3 


-5 


-4.4 




NGC 4889 


E4 


-5 


-4.1 




NGC 3245 


SO 


-2 


-4.6 




NGC 3941 


SBO/a 





-2.7 




NGC 4262 


SBO 


-2 


-3.8 




NGC 5866 


SO 


-2 


-4.4 




NGC 1357 


Sa 


1 


-2.0 




NGC 2775 


Sa 


1 


-3.4 




NGC 3368 


Sab 


2 


-3.3 


L2 


NGC 3623 


Sa 


1 


-4.1 


L2: 


NGC 1832 


SBb 


3 


2.4 




NGC 3147 


Sb 


3 


0.0 


S2 


NGC 3627 


Sb 


3 


1.4 


T2/S2 


NGC 4775 


Sc 


5 


9.7 




NGC 5248 


Sbc 


4 


1.4 


H 


NGC 6217 


SBbc 


4 


4.8 


H 


NGC 2903 


Sc 


5 


1.1 


H 


NGC 4631 


Sc 


5 


9.4 


H 


NGC 6181 


Sc 


5 


3.1 


H 


NGC 6643 


Sc 


5 


3.8 


H 


NGC 4449 


Sm/Im 


9 


10.8 


H 



for the values of the principal components. However, it 
can be shown that, to a large extent, these solutions are 
equivalent, all producing essentially the same ranking of 
galaxies along the first principal component. 

Fig. 1 shows the projection of the spectra of the 23 
standards on to the principal plane. This plane is highly 
informative about the structure of the spectral space, as 
it contains 94.8% of the data variance. The spectral se- 
quence is readily identified along the first principal com- 
ponent. We then define quantitatively the spectral type 
ST of a galaxy by the value of its first principal compo- 
nent. These values are presented in table 1. Note that ST 
increases from early to late types. In what follows, the 
spectral sequence will refer to that obtained with the 23 
standards of table 1. The mean values and dispersions of 
ST for the standards are -4.9 ± 0.8 for E-E/S0, -3.9 ± 0.9 
for SO-SO/a, -3.2 ± 0.9 for Sa-Sab, 2.0 ± 1.8 for Sb-Sbc, 
and 6.3 ±4.1 for Sc-Im. The ranking of spectral types in 
these groups is analogous to the ranking of Hubble mor- 
phological types. Note that the dispersion of ST around 
the mean value of each group is large enough to produce 
a substantial overlap in the spectral sequence of galaxies 
of different Hubble types: galaxies of same morphologies 
may have very different spectral types. 

It is worth emphasizing that we base our spectral clas- 
sification only on the integrated properties of the stellar 
populations that are contained in the continuum and ab- 
sorption lines, and that the emission lines enter in no way 
in the classification scheme. 
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Fig. 1. Projection of the spectra of the 23 standard galax- 
ies on to the principal plane. The spectral sequence follows 
the first principal component, and a spectral type is at- 
tributed to each galaxy by its value of this component. 
Filled symbols: E and SO galaxies; open symbols: S galax- 
ies. 

2.3. Emission line data and reddening correction 

Among our 23 standard galaxies, emission lines are mea- 
sured only for the 15 spirals (T > 1), being too weak or 
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undetected in the integrated spectra of galaxies of ear- 
lier types (at the resolution of the observations, ~5-7A; 
K92b). We use all the emission lines with equiva- 



see 



lent widths (EW) measured by jK92a} [OII]A3727, H/3, 



[OIILU5007, Ha, [NII]AA6548, 6583, and [SII]AA6717, 
6731. The EW of Ha was computed from the EW of 



Ha+[NII] and the value of [Nil] /Ha, also from |K92a . 
Other emission lines, like [OI]AA6300,6363, are not strong 
enough to b e relia bly measured in the spectra of Kenni- 
cutt's atlas ( |K92bD . 

In the next sections we will discuss the behavior of red- 
dening corrected line ratios. For that, we have obtained 
a value of the effective internal reddening of the galax- 
ies from an estimation of the emission line intensity ratio 
of Ha and H/3, Ha/H /3. Our procedure is similar to the 
one adopted by K92a . We begin by modeling the relation 
between the observed EWs of Ha and H/3 by a straight 
line (Fig. 2). Usin g the bisector estimation for the fitting 
(Isobe et al. 1990), we have obtained 



EW(H/3) ofcs = -4.19(±0.36) + 

0.194(±0.011)EW(Ha) ofcs 



(1) 



Note that this linear relationship is purely empirical, and 
there is not, a priori, any reason for EW(H/3) h s and EW 
(Ha) h s be linearly related. 

The observed EW of a line may be written as EW i s 
= EW em + EW^, where EW em is the emission compo- 
nent (> 0) and EW a ;, s corresponds to the component due 
to stellar absorption (< 0). To proceed, we need to know 
how EW a i, s behaves as a function of EW em or EW f, s . Let 
us suppose, then, that the emission and observed EWs of 
Ha and H/3 may also be related linearly, what is a rea- 
sonable assumption in most models (e.g., Barbaro & Pog- 
gianti |1997| ). Writing EW(Ha) em = a a + 6 Q EW(Ha) obs 
and EW(H/3)em — ap + 6 ( 3EW(H/3) fc s , and substituting 
these relations in Eq. 1, we obtain EW(H/3) em = C + 
R EW(Ha) em , where C and R are constants. Now, physi- 
cal solutions require C equal to zero, because EW(H/3) em 
should go to zero when EW(Ha) em goes to zero. Conse- 
quently, we obtain that the emission EWs of Ha and H/3 
are related by 



EW(H/3) em = R EW(Ha) e 



(2) 



with R = 0.194& ( 3/& a . With these assumptions, Ha/H/3 
depends only of R and the ratio of the continuum fluxes 
at Ha and H/3, that we have measured in the spectra. 
We have adopted here R = 0.2, since models like those of 
Barbaro & Poggianti (1997) indicate that bp ~ b a . 

The extinction coefficient at H/3, C(H/3), has been com- 
puted assuming an intrinsic value for Ha/H/3 of 2.9 (cor- 
responding to an electron temperature of 9000K in the 
case B of the theory of recombination-line radiation). The 
expected formal error in C(H/3) from this procedure is 
~ 0.18. The values of C(H/3) may be considered a crude 
estimation of the global extinction in the galaxies, since, 



21 




EW(Ha) (A) 

Fig. 2. Relation between the observed values of the equiv- 
alent widths of Ha and H/3. The straight line corresponds 
to the fitting given by Eq. 1. 

as mentioned above, foreground extinction by our Galaxy 
is negligible. 

It is important to stress that, although very simple, 
this estimate of the reddening is quite robust, and there is 
not any evidence in the data that a more complex model 
is required. As we shall see, it leads to a Ha/H/3 ratio 
that decreases steadly with spectral type. We have verified 
that this trend continues when other ways to estimate the 
reddening are used. Also, the behavior of other line ratios 
discussed in the next sections are, to a large extent, not 
strongly affected by the use of other values of R. 

3. Trends of emission line properties with galaxy 
spectral type 

3.1. Equivalent widths 

Fig. 3 consists of 5 panels showing the behavior of the ob- 
served EWs of [Oil], H/3, [OIII], Ha+[NII] and [SII] as a 
function of ST. In all of them we see a clear increase of 
the emission line EWs as a function of the galaxy spectral 
type. For [OIII], the increase is apparent only from ST=2 
upwards, since this line is too weak to be effectively mea- 
sured in spectral types lower (i.e., earlier) than that. We 
had shown, in Fig. 2, the relation between EW(H/3) and 
EW(Hq). Fig. 4 shows the relation of the EWs of [Oil], 
[OIII], and [SII] with EW(Ha+[NII]). The EWs of all 
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Table 2. Extinction and corrected line intensities of spiral galaxies 



Name* 


C(H/3) 


[OII]/H/J 


[OIII]/H/3 


[Nil] /Ha 


[SII]/Ha 


Ha/H/3 


NGC1357 


0.71 


1.70 


<0.14 


0.23 


0.05 


5.12 


NGC2775 


0.72 


3.10 


<1.11 


0.28 


- 


5.20 


NGC3368 


0.73 


4.76 


<0.74 


0.47 


0.16 


5.22 


NGC3623 


0.72 


3.93 


<1.79 


0.83 


0.27 


5.18 


NGC1832 


0.52 


2.21 


0.26 


0.15 


0.07 


4.42 


NGC3147 


0.63 


0.94 


<0.18 


0.16 


- 


4.83 


NGC3627 


0.58 


1.05 


<0.08 


0.16 


0.06 


4.65 


NGC4775 


0.16 


3.40 


1.25 


0.08 


0.10 


3.31 


NGC5248 


0.58 


1.28 


0.13 


0.14 


0.06 


4.62 


NGC6217 


0.40 


1.49 


0.26 


0.18 


0.07 


3.99 


NGC2903 


0.64 


1.29 


<0.09 


0.16 


0.05 


4.84 


NGC4631 


0.24 


3.44 


1.20 


0.07 


0.10 


3.52 


NGC6181 


0.55 


1.51 


0.30 


0.14 


0.06 


4.51 


NGC6643 


0.49 


1.35 


0.18 


0.13 


0.07 


4.31 


NGC4449 


0.16 


3.52 


1.82 


0.04 


0.08 


3.31 


a The sign 


:' after the 


galaxy name 


indicates that EW(H/3) 


obs < 





these lines correlate well with EW(Ha) or EW(Ha+[NII]). 
The EWs of the Balmer emission lines are considered good 
tracers of the star formation rate in galaxies (Kennicutt 



1983), and it has been shown that [Oil] is a good sub- 



stitute of the Balmer lines in high-redshift galaxies (Gal- 
lagher et al. 1989; K92a). Our results indicate that, for 



practical purposes, any of the EWs displayed in Fig. 4 (ex- 
cept maybe [OIII] , which is available for a smaller range of 
galaxy types) can be used to estimate the star formation 
rate in normal galaxies, after empirical calibration on the 
Balmer line EWs. It is worth noting, however, that deriv- 
ing a star formation rate from an emission line equivalent 
width is model dependent because, among other things, 
such a procedure assumes that the ionization is provided 
by the radiation field from massive stars, which is not nec- 
essarily true, specially in early type galaxies (see below). 
On the other side, the EW of any of these lines is an em- 
pirical indicator of galaxy spectral types, given the good 
correlation between these two quantities, shown in Fig. 3. 
It is interesting to compare our Figs. 2 and 4 with Figs. 



4 to 6 of K92a, which include all galaxies with emission 



lines in K92b. Most of the scatter observed in the K92a 



plots is produced by galaxies not belonging to the sample 
of normal galaxies discussed here. 

3.2. Emission line intensity ratios 

Fig. 5 consists of 9 panels showing the behavior of several 
emission line ratios as a function of the galaxy spectral 
type. Line ratios are corrected for reddening, as explained 
in Sect. 2.3. 

Fig. 5 shows that very clear trends exist between emis- 
sion line ratios and spectral types, with very small scatter. 
Briefly, the trends can be described as follows. 

1. Ha/H/3 (panel 5-a) decreases steadily with increasing 
spectral type. Converted into total extinction at H/3, 



this means that C(H/3) decreases from about 0.7 at ST 
—5 to 0.2 at ST - +10. 

2. [OII]/H/3 (panel 5-b) first decreases with spectral type, 
for ST ^ -2, then increases slowly. 

3. [OIII] is measurable only from ST ~ 1 upwards, where 
[OIIIJ/H/3 (panel 5-c) increases very strongly with ST. 
For ST 1 only upper values are available for the 
[OIII] emission but [OIII]/H/3 seems to decrease in the 
interval -5 <, ST < 0. 

4. [Nil] /Ha (panel 5-d) decreases abruptly for ST < 0, 
decreasing mildly towards larger values of ST. 

5. [SII]/Ha (panel 5-c) also decreases for ST < 0, then 
increases smoothly for larger values of ST. 

6. ( [Oil] + [OIII] )/H/3 (panel 5-f) first decreases until ST 
~ (this effect is due to [Oil] only since [OIII] is weak 
in this range of ST) , then increases steadily, but with a 
smaller slope than [OIII]/H/3. Note that, contrarily to 
what one might think, the values of ([OH] + [OIII]) for 
early-type galaxies are not upper values, since in this 
case the [Oil] emission is much larger than the [OIII] 
emission. 

7. [OIII] /[Oil] (panel 5-g) increases steadily for ST > 1 
(the effect of [OIII] is dominant). For lower values of 
ST only upper values of [OIII] /[Oil] are known, but 
there is a hint that this ratio decreases with ST. 

8. [Nil] /[Oil] (panel 5-h) is roughly constant or decreases 
slowly until ST ~ 5, decreasing faster for larger values 
of ST. 

9. [SII]/[OII] (panel 5-i) shows little trend with respect 
to ST, except, maybe, a slight decrease with increasing 
spectral types. 

The trends seen in Figs. 3 and 5 are impressive, con- 
sidering that the galaxies were selected simply using a 
criterion of "normality", and that the ordering by spec- 
tral type relies only on the continuum and absorption fea- 
tures of the integrated spectra. However, if we recall that 
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lines arise due mainly to the presence of hot stars (T e ff 
30000K) in the stellar populations, and are modulated 
by other parameters such as chemical abundances and 
gas content, which also seem to be related to the gen- 



eral galaxy properties (Roberts & Haynes 1994). Individ- 
ual HII regions in a galaxy show a definite scatter in their 
properties (e.g., Kennicutt & Garnett 1996| , Roy & Walsh 
p.997 ), but their total number in a galaxy is very large 
(typically of the order of one hundred) so that the inte- 
grated emission line spectrum of a galaxy is well defined. 
The fact that we find a close link between the emission 
lines and the continuum and absorption features that are 
behind the spectral classification suggests that there is a 
simple relation between the population of hot, ionizing 
stars and the population of stars responsible for the opti- 
cal continuum. 

It is interesting to see how the emission line ratios be- 
have as a function of the de Vaucouleurs type of the galaxy, 
which is so far the most widespread way of ordering galax- 
ies and is based on their morphological type. Fig. 6 shows 
the same line ratios as Fig. 5, but as a function of the 
de Vaucouleurs T type. One can distinguish trends simi- 
lar to those seen in Fig. 5, but they are much less clear. 
That similar trends can be seen is not surprising, since 
there is some correlation between the spectral sequence 



and the Hubble sequence (Sodre & Cuevas 1997, Folkes 
et al. 1996). That the trends are less clear than those in 
Fig. 5 can be explained by assuming that the populations 
of hot stars are not expected to be correlated with galaxy 
morphological types to the same extent as they are ex- 
pected to be correlated with the entire stellar populations, 
as measured by the spectral types. 



in 
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Fig. 4. Relation of the observed equivalent widths of a) 
[Oil], b) [OIII], and c) [SII] with EW(Ha+[NII]). 



these features reflect the galactic stellar populations, and 
if the stellar populations obey some simple rules, then it 
is not so surprising that the emission line properties from 
integrated spectra can also be ordered. Indeed, emission 



4. Discussion 

Before trying to interpret the trends shown above, it is 
useful to recall some basics on emission lines and pho- 
toionization models. This is done in the next subsection. 



4-1- A reminder on emission line theory 

Schematically, the (reddening corrected) intensity ratios 
of emission lines produced in photoionized regions are a 
function of the following parameters (see e.g. Stasihska 



1998 and references therein): 



1. the global metallicity, O/H; 

2. abundance ratios, like N/O and S/O; 

3. the mean effective temperature of the ionizing radia- 
tion field < T e f f > ; 

4. the average ionization parameter U = Qh /(^R^n), 
where Qh is the number of ionizing photons of an HII 
region, R s is the Stromgren radius, and n is the gas 
density. 

The effect of these parameters on the line ratios are 
the following. 
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Oxygen is the major coolant in HII regions. As O/H 
increases, the electron temperature decreases because 
of the increased cooling of the gas. The [OIIIJ/H/3 ra- 
tio first increases (due to an abundance effect), then, 
for O/H greater than about half solar, decreases due 
to an electron temperature effect: the gas becomes so 
cool that the optical [Offf]A5007 line gets difficult to 
excite. In this case, cooling occurs through the [OIII] 
far infrared lines at 52/xm and 88/im. Qualitatively, the 
[OII]/H/3 ratio behaves similarly to [OIII]/H/3, but it is 
not so reduced at high metallicities, because the cool- 
ing is less efficient in the region emitting [OH] than 
in the region emitting [OIII]. Line ratios like [Nil] /Ha 
or [SII]/Ha are affected by O/H through the electron 
temperature: they get enhanced as O/H decreases. 
An increase of N/O only enhances [Nil] /[Oil], but does 
not affect the other optical line ratios, since nitrogen 
contributes little to the cooling at the abundances ex- 
pected for the general interstellar medium. Similarly, 
an increase of S/O enhances [SII]/[OII] and leaves the 
rest unchanged. However, one does not expected S/O 
to vary among galaxies, since sulfur and oxygen are 
produced in the same stars. Therefore, any change in 
[SII]/[OII] should rather be attributed to other causes. 
The effects of < T e // > are twofold. Firstly, < T e ff > 
acts on the ionization structure: the proportions of + , 
N + or S + ions decrease with increasing < T e ff >. Sec- 
ondly, it influences the thermal balance of the nebula: 
as < T e ff > increases, the energy gains become larger 
and the electron temperature rises, increasing the in- 
tensities of the forbidden lines with respect to Ha or 
H/3. 

The effects of decreasing U are to reduce the average 
ionization, and to decrease ratios like ++ /0 + . For 
very low U , such as found in the diffuse interstellar 
medium, a wide transition zone of low ionization de- 
velops, containing 0° and S + ions and still hot enough 
to allow collisional excitation of optical forbidden lines. 
Therefore, the [SII|/Ha and the [SII]/[OII] ratios are 
higher than in HII regions having a large U. 
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Fig. 7. Classical emission line ratio diagnostic diagrams. 
The data for galaxies of our sample (circles) are plotted 
together with data f or the giant HII regions observed by 
McCall et al. fll985| ) (dots). Upper values (of [OIII]) are 
indicated by arrows. The panels are: a) [OIII]/H/? ver- 
sus [OII]/H/3, b) [OIII] /H/3 versus [NII]/Ha, c) [Nil] /Ha 
versus [OII]/H/?, d) [OIII]/H/3 versus [SII]/Ha, and e) 
[SII]/Ha versus [OIIJ/H/3. 



If shocks are present, the emission line ratios are mod- 
ified. The effect of shocks is to heat the gas to very high 
temperatures (T e = I0 6 — 10 7 K). By recombination and 
free-free processes, this produces a hard radiation field 
which strongly heats the post shock region, producing an 
extended, warm, transition region, and enhancing the lines 
that are most sensitive to the temperature. As a result, 
[OII]/H/3, [Nil] /Ha, and [SII]/Ha line ratios are enhanced 
with respect to pure photoionization nebulae. 

4-2. Diagnostic diagrams 



Line rati o diag rams (Baldwin et al. 1981, Veilleux & Os- 
terbrock 1987 ) are helpful for the diagnostics of emis- 
sion line objects. We now use such diagrams to compare 



the emission properties of galaxy integrated spectra with 
those of individual giant HII regions. Figs. 7 and 8 show 
the data from the integrated spectra of the normal spi- 
ral galaxies together with those of the giant HII regions 
(GHRs) observed by McCall, Rybski & Shields flL985D in 
several spiral galaxies and which compose the GHR se- 
quence. This sequence is interpreted as being a sequence in 
mctallicity (O/H) and also in mean effective temperature 
or in mean ionization parameter (or both) (McCall et al. 
|f 985 , Dopita & Evans 198E). Note that radial abundance 
gradients, which are a common feature in spiral galaxies 



(Zaritsky et al. 1994), complicate the interpretation of in- 
tegrated spectra: the contribution of each annular region 
is weighted by the luminosities of the HII regions found 
there and by their spectral properties. 
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Panel 7-a shows [OIII] /H/3 versus [Oil] /H/3. We see 
that the points corresponding to the integrated spectra of 
our standard galaxies that have [OIII] large enough to be 
observed (i.e. those with ST >0) are well inside the GHR 
sequence. 

Panel 7-b displays [OIII] /H/3 versus [Nil] /Ha. Galax- 
ies with measured [OIII] emission are disposed along the 
GHR sequence. The outliers are galaxies for which only 
upper values of the [OIII] emission are available. They are 
early type spirals, with ST < 0, as can be deduced from 
Figs. 5c and 5d. These galaxies stand out conspicuously in 
the [Nil] /Ha versus [Oil] /H/3, diagram (Panel 7c) as well. 
As will be seen below, we interpret this behavior as being 
the combined effect of an overabundance of nitrogen in 
early-type spirals and an ionization source different from 
ordinary O stars. 

Panel 7-d displays [OIII]/H/3 versus [SII]/Ha. In this 
diagram, while the standard galaxies with measured log( 
[OIII]/ H/3) < 1 are within the HII region sequence, those 
with log( [OIII]/ H/3) >1 (i.e. all the galaxies with ST > 
9) are slightly to the upper right. This may indicate an in- 
creasing contribution of a diffuse ionized medium with in- 
creasing galaxy spectral type. Note that Lehnert & Heck- 



man (1994|), plotting the spectra of K82a having EW(Ha) 
> 30 Ah [OIII]/ H/3 versus [Nil]/ Ha and [OIII]/ H/3 
versus [SII] / Ha planes also concluded that diffuse ionized 
gas contributes to the integrated spectra of galaxies. [] 

Panel 7-e shows [Oil] /H/3 versus [SII] /Ha. The same 
early-type spirals that show enhanced [Nil] emission also 
show enhanced [SII] e mission compared with the GHRs 
sample of McCall et al. 1985 , although to a lesser extent. 
They are among the ones with the largest [Oil] /H/3 ratios. 
We will argue below that the bulk of the emission in these 
objects is not due to ordinary O stars. 

Fig. 8 shows some line ratios as a function of ([OII] + 
[OIII])/ H/3. Since the pioneering work of Pagel et al. 
( |l979| ), the ([OII]+ [OIII]) /H/3 ratio has been widely used 
to derive the oxygen abundance. In principle, the relation 
between ([ Oil] + [0111])/ H /3 and O/H is double sided (e.g. 
McGaugh 1991, Stasihska 1998] ). In integrated spectra of 
galaxies, regions within a few kiloparsecs from the galactic 
center have probably the greater weigh t, bec ause most of 
the star formation is there (Rozas et al. 1996 ). Since these 
regions are the most metal rich, one expects that the in- 
tegrated spectra correspond to the regime where ([OII] + 



1 The detection of a diffuse medium is actually not pos- 
sible using a [OIII]/ H/3 vs [ NH]/ Ha diagram. Indeed, pho- 



toionization models (Stasiriska 1998) shows that decreasing the 
ionization parameter alone to very low values does not shift 
the points out of the GHR sequence. On the other hand, this 
diagram is sensitive to the N/O ratio which, unlike S/O, is 
expected to vary significantly from object to object. Detailed 
spectroscopic observations of t he di ffuse intergalactic medium 
in spiral galaxies (Wang et al. 1997 ) however suggest that the 
ionization of this medium cannot be entirely attributed to ra- 
diation from O stars. 



[OIII])/ H/3 increases as O/H decreases (but this should 
be confirmed by simulations). 

Panel 8-a shows [OIII] /[Oil] versus ([OII]+ [OIII])/ 
H/3. [OIII]/ [Oil] roughly indicates the ionization state, 
which is linked to the ionization parameter U, thus to the 
gas density distribution in each HII region, and to the 
hardness of the ionizing radiation field. In this diagram, 
the data points for the integrated galaxy spectra lie well 
within the region occupied by GHRs except, maybe, for 
the early spectral type galaxies, with only upper values of 
[OIII]/ H/3, that tend to locate in the lower envelope of 
the GHR sequence. The spread of ([OII]+ [OIII])/ H/3 for 
the galaxies is smaller than for individual GHRs. This is 
reasonable, since the integrated spectra are averages over 
GHRs of various metallicities, ionization parameters and 
mean affective temperatures. 

Panel 8-b displays [Nil] /[OR] versus ([OII]+ [OIII])/ 
H/3. We see that all galaxies in the sample fall inside the 
GHR sequence in fig. 8-b, except three that have higher 
[Nil]/ [Oil] than GHRs of same ([OII]+ [OIII])/ H/3. Note 
that the three exceptions also have high [Nil] /Ha and are 
all the galaxies of our sample that have ST < -3. 

Panel 8-c displays [SII] /[Oil] versus 
([Oil] + [0111]) /H/3. In this diagram, most of the standard 
galaxies fall inside the GHR sequence (although, as 
already noted by McCall et al. 1985, the GHR sequence 



is not so well defined when using [SII] /[Oil] instead of 
[Nil] /[Oil]), the galaxies with ST < -3 tending to be 
above that sequence. 

4-3. A tentative interpretation of emission line trends 
with galaxy spectral types 

While the trends between the emission line properties in 
the integrated spectra of normal spiral galaxies and spec- 
tral types are impressive, their interpretation is not neces- 
sarily straightforward, since it is dependent of the physical 
state and spatial distribution of the gas and of the ioniza- 
tion mechanisms. 



It was shown by Zaritsky et al. (1994) that the char- 
acteristic oxygen abundance in a spiral galaxy derived 
from its individual giant HII regions increases towards 
earlier morphological types. The trend seen in Fig. 5-f 
for ST > corresponds (at least qualitatively) to what 



is expected. Using the Zaritsky et al. (1994) calibration 



of ([Oil] + [0111]) /H/3 into O/H, it would translate into a 
decrease in the average O/H by ^0.5 dex for ST going 
from to 10. This is only an indicative value, because the 
presence of abundance gradients makes it impossible to 
be more informative without numerical simulations. The 
case of the galaxies with ST < 0, will be discussed later. 



Zaritsky et al. (1994) have also found a strong corre- 



lation between the characteristic oxygen abundance and 



galaxy luminosity ( see al so Vila-Costas & Edmunds 1992 



Roberts & Haynes 1994 ). Due to the strong correlation 
between galaxy luminosity and Hubble type in their sam- 
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pie, it is impossible to know which of these two quantities 
is the primary cause of the correlation. It is interesting to 
point out that, in the sample of galaxies discussed here, 
there is no such correlation between O/H and luminosity 
(most of the galaxies in table 2 have in the narrow 
range between -19.9 and -20.7, see K92b| ) and that the cor- 
relation between luminosity and spectral type, if existent, 
is at most weak. This suggests that O/H is rather linked to 
the stellar populations, as measured by the spectral types, 
than to galaxy luminosity. 

The increase of [OIII]/[OII] with spectral type (Fig. 5- 
g) for ST ^ is probably due to an increase of < T e f f > 
(or an increase of U) when O/H decreases, similarly to 
what is invoked for individual GHRs (e.g., Garnett & 



Shields |1987|) . 

The increase of [Nil] /[Oil] with decreasing spectral 
type (Fig. 5-h) is probably due to an increase of N/O. 
Indeed, in GHRs, the fact that [NII]/[OII] increases with 
decreasing ([Oil] + [0111]) /H/3 is interpreted, using pho- 
toionization models, as being due to an increase of N/O 
with O/H (Thurston, Edmunds & Henry |1996| ). It is likely 
then, that in the spectral sequence of galaxies there is 
an increase of the average N/O with decreasing spectral 
type, linked with the increase of O/H. This is in qualita- 
tive agreement with a secondary production of nitrogen. 
A priori, one cannot exclude the fact that the increase of 
[Nil] /[Oil] with decreasing spectral type could be simply 
attributed to the decrease of O/H (i.e., with N/O stay- 
ing constant), due to the thermal properties of GHRs. In- 
deed, [Nil] /[Oil] = (N/O) (e(NII)/e(OII)), where e(NII) 
and e(OII) are the emissivities of the [Nil] and [Oil] lines, 
and e(NII)/e(OII) is a decreasing function of the electron 
temperature, thus an increasing function of O/H. Note, 
however, that the three early-type galaxies that are out of 
the GHR sequence in Figs. 7-c, 7-e, and 8-b, must have 
particularly high N/O, because their large [Nil] /[Oil] can- 
not be due to a particularly low electron temperature, 
since [Nil] /Ha is large as well. 

The ratio [SII]/[OII] does not present any significant 
trend with spectral type (Fig. 5-i). This is consistent with 
the assumption that the abundance ratio S/O is constant, 
since both S and O are primary elements. There is a hint of 
a decrease in that line ratio with increasing spectral type 
but this may be ascribed to the variations of the ionization 
parameter along the sequence more than to variations in 
the relative abundances of S and O. 

The increase of [Oil] /H/3 as the galaxy spectral type 
decreases from downwards (as well as that of [SII]/Ha 
and [Nil] /Ha, Figs. 5-b, d, e) is very interesting. Early 
spectral type galaxies appear as upper limits in the [OIII]/ 
H/3 versus [Oil] /H/3 diagnostic diagram, since their [OIII] 
emission is weak. They stand out completely from the gi- 
ant HII region sequence in diagnostic diagrams with [Nil] 
or [SII] . The integrated spectra of these galaxies are sim- 
ilar to those of LINERs. 



Interestingly, 11 of the 15 galaxies of our sample 
have had their nuclear regions observed by Ho et al. 
( 1997| ). Their classification of these nuclei into LINERs 
(L), Seyfert 2 (S2) and HII regions (H) is as follows (see 
table 1): all the H type nuclei correspond to ST > 0, while 
NGC3623 and NGC3368 (ST < -3) correspond to LIN- 
ERs. NGC3147 (ST = 0.0) has a S2 nuclear spectrum and 
NGC3627 (ST = 1.4) has a transition spectrum (T2/S2). 
The other gala xies in table 1 have not been observed by 
Ho et al. ( 1997 ). So, the integrated spectra of the galaxies 
with earliest spectral types have characteristics of LIN- 
ERs, and the nuclear regions of these galaxies too (we have 
checked that the contribution of the LINER nuclei to the 
integrated galaxy spectra is negligible). This seems to in- 
dicate that what gives rise to the LINER phenomenon, at 
least in the normal galaxies we are studying here, is not 
specifically related to the nucleus. 

Given the evidence presented above that metallicity, as 
measured by O/H, decreases along the spectral sequence, 
it would be tempting to attribute the distinct behavior of 
the galaxies having ST < in Figs. 5, 7, and 8 to over- 
abundances in the gas. But, as stated in Sect 4.1, over- 
abundances would produce low [Oil] /H/3, not large ones 
as observed. 

Another possibility is that the distinctive features pre- 
sented by early type galaxies are due to the aging of the 
stellar populations associated with GHRs. Sodre & Cuevas 



(1997), using the spectral synthesis code GISSEL (Bruzual 
& Chariot 1995), have shown that the spectral sequence 
of galaxies is well reproduced if one assumes that the star 
formation rate has the form exp(— i/£*), where is an in- 
creasing function of the galaxy spectral type. Such a model 
explains, at least qualitatively, why the equivalent width 
of H/3 emission increases with increasing galaxy spectral 
type, as found in Sect. 3. For galaxies with the earliest 
spectral types, most of the star formation has occurred 
long ago, and very few O stars are present. In this case, 
the radiation field in the Lyman continuum may be dom- 



inated by post-AGB stars (e.g. Bressan et al. 1994), be 



cause, in simple starburst models, these stars provide most 
of the ionizing photons for bursts older than ~ 10 s yr. This 
radiation field is much harder than that produced by the 
massive O stars that power classical HII regions. It is avail- 
able to ionize the gas bound in the primitive HII regions 
or the diffuse material produced by their eventual disrup- 
tion. Binette et al. ( |1994 ) have shown that the emission 
line properties (EWs and line ratios) of early type galaxies 
can be accounted for by photoionization by such a stellar 
population. A similar model could well explain at the same 
time the high [Oil] /H/3, [SII] /Ha, and [Nil] /Ha found in 
the early spectral type galaxies of our sample. 

The long lasting dispute about the ionization mech- 
anism in LINERs - photoionization versus shocks - is 
not settled yet, precisely because under certain conditions 
both explanations can account for the observed emission 
line spectra. Models of evolving stellar populations that es- 
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timate the energy release in supernovae and stellar winds 
(e.g., Leitherer & Heckman 1995| ) will help in the solution 
of the problem. 

We have seen that the effective extinction coefficient 
at H/3, C(H/3), decreases steadily with increasing spectral 
type. The beautiful trend shown in Fig. 5-a is somewhat 
surprising, since it is not commonly believed that early- 
type spirals have higher average extinction than galaxies 
of later types. Indeed, the related issue of the overall opac- 
ity in galaxies has been subject of significant debat e, often 
with contra dictor y opinions (e.g., Valen tjin 199C , Davies 
& Burstein 1995 ). Wang & Heckman ( 1996 ) have veri- 
fied that the dust opacity, measured from the ultraviolet 
to far-infrared luminosity ratio, increases with increasing 
galaxy infrared luminosity. Taking into account that the 
far-infrared luminosity is approximately constant for early 
morphological type spirals and decreases slowly towards 
later morphological types (e.g., Roberts & Heynes 1994 ), 
the correlation noticed by Wang & Heckman might indi- 
cate that the dust opacity is indeed decreasing towards 
later morphological types, in qualitative agreement with 
the results reported here. 

Now, assuming for simplicity that dust properties are 
invariant with galaxy spectral type, this trend suggests 
that the mean surface density of dust in galaxies decreases 
along the spectral sequence from early to late spectral 
types by a factor of 4. Since the mean surface density 
of neutral hydrogen in normal galaxies increases towards 
later types (Roberts & Heines 1994), the dust-to-gas ratio 



must be decreasing as the spectral type increases. This 
is consistent with the idea of dust being formed less effi- 
ciently in low metallicity environc mcnts , as shown obser- 
vationally by Lisenfeld & Ferrara ( 1998 ). Note that some 
of the scatter in Fig. 5-a may be produced by the random 
inclination of the galaxies and that the far infrared prop- 
erties of normal galaxies actually suggest some variations 
of th e dust properties with galaxy type (Sauvage & Thuan 



1994|) 



5. Summary 

We have analyzed the global emission line properties of 
normal spiral galaxies from their integrated spectra. Ini- 
tially, we have selected from the Kennicutt's atlas (K92b) 
23 galaxies which do not present any sign of spectral or 
morphological peculiarity. A PCA based only on continu- 
ous and absorption features of the integrated spectra of 
these galaxies enabled us to define the galaxy spectral 
types. We have then considered the subsample of 15 galax- 



ies whose integrated spectra in K92b revealed emission 



lines. All these galaxies are spirals. We have studied the 
progression of the properties of these emission lines with 
galaxy spectral type. 

Our major results are the following: 

1. The EWs of all the lines correlate well with EW(Hct). 
As a consequence, for practical purposes, the EW of 



any of the studied lines ([OII]A3727, H/3, [OIII]A5007, 
Ha, [NII]AA6548, 6583, and [SII]AA6717, 6731) can be 
used to estimate the average massive star formation 
rate, at least for galaxies with spectral types <; 0. 
Galaxies with lower spectral type have emission line 
spectra resembling low excitation LINERs and, conse- 
quently, estimates of their global star formation rate 
based on the their emission lines are unreliable. 

2. The EW of all lines correlate very well with ST, sug- 
gesting that they can be used as indicators of galaxy 
spectral types. 

3. Emission line ratios present a remarkable regularity 
when plotted against galaxy spectral types. This stems 
from the fact that the strengths of the emission lines 
are linked to the stellar populations, and the latter con- 
stitute the basis of the spectral classification of galax- 
ies. The plots are much noisier when using ordinary 
morphological types. 

4. For galaxies with ST > 0, the increase of (OII+ OIII)/ 
H/3 with increasing spectral type suggests a decreasing 
average metallicity along the spectral sequence. 

5. The increase of [OIII] /[Oil] with increasing spectral 
type may be due either to an increase of the mean ef- 
fective temperature of the ionizing stars, < T e ff >, or 
to an increase of the characteristic ionization parame- 
ter U. 

6. The increase of [Nil] /[Oil] with decreasing spectral 
type is probably due to an increase of the average N/O. 
This behavior is consistent with the fact that, in giant 
HII regions, N/O increases with O/H. 

7. The ratio [SII]/[OII] is essentially constant with spec- 
tral type and is consistent with S/O being constant. 
This is expected because S and O are both primary 
elements. 

8. The integrated spectra of early-type spirals (ST ^ 1) 
are similar to those of LINERs. The long character- 
istic time for star formation in e arly-type spirals as 
derived by Sodre & Cuevas ( 1997 ) provides a natural 
explanation in terms of ionization by post-AGB stars. 
Interestingly, in several of these galaxies, the nuclei 
themselves have been ch aracte rized as dwarf LINERs 
in the work of Ho et al. ( 1997 ). 

9. The extinction by dust in HII regions seems to decrease 
with increasing spectral type. This may be a direct 
consequence of the decrease in mean metallicity along 
the spectral sequence. 

While the interpretation of the trends found with spec- 
tral type is only tentative, the reality of these trends is 
beyond any doubt. 

This work confirms that the spectral classification is 
useful in the analysis of galaxy properties, being able to 
replace, with advantages, the traditional morphological 
classification for many purposes. It produces an objec- 
tive, quantitative and easy to implement classification, 
that may be more satisfactory for the analysis of large 
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data bases than the subjective and hard to obtain mor- 
phological types. Also, since spectra are closely linked to 
the physics of the stars and of the interstellar medium 
that build up galaxies, spectral types are more adequate 
than morphological types for a quantitative approach in 
studies of galaxy properties, as shown by this work. 
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Fig. 3. Relation between the observed equivalent widths of a) [Oil], b) H/3, c) [OIII], d) Hq;+[NII], and e) [SII] and 
galaxy spectral type (ST). Upper values of EW([OIII]) are indicated by arrows in panel c. 
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Fig. 5. Behavior of emission line ratios measured in the galaxy integrated spectra as a function of the galaxy spectral 
type (ST). The line ratios (except Ha/H/J) are corrected for reddening as described in the text. The emission line ratios 
in the Fig. are: a) Ha/H/3, b) [OII]/H/3, c) [OIIIJ/H/3, d) [Nil] /Ha, e) [SII]/Ha, f) ([OII] + [OIII])/H/3, g) [OIII]/[OII], 
h) [Nil] /[Oil], i) [SII]/[OII]. Arrows indicate upper values. 
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Fig. 6. Same as Fig. 5, but now the line ratios are plotted as a function of the de Vaucouleurs type of the galaxies. 
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Fig. 8. Forbidden line ratios as a function of the O/H indicator ([OII] + [OIII])/H/3. The dat a for galaxies of our sample 
(circles) are plotted together with data for the giant HII regions observed by McCall et al. ( 1985 ) (dots). Upper values 
(of [OIII]) are indicated by arrows. The panels are: a) [OIII]/[OII] versus ([OII] + [OIII])/H/3, b) [Nil] /[Oil] versus 
([OII] + [OIII])/H/3, and c) [SII]/[OII] versus ( [Oil] + [OIII] )/H/3. 



